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In this note we present a scaling theory for the elasticity of Olympic gels, i.e., gels where the 
elasticity is a consequence of topology only. It is shown that two deformation regimes exist. The 
first is the non affine deformation regime where the free energy scales linear with the deformation. 
In the large (affine) deformation regime the free energy is shown to scale as F oc A®^'^ where A is the 
deformation ratio. Thus a highly non Hookian stress - strain relation is predicted. 
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In this communication we compute the scaling of the 
elastic modulus of an Olympic gel. This name has been 
created by de Gennes [|] since the configuration of such 
gels resemble very much the Olympic rings. Olympic 
gels are very peculiar elastic materials as their elasticity 
does not come from crosslinks as in conventional rub- 
bers. Actually ideal Olympic gels do not contain any 
crosslinks, but consist of concatenated rings only. In this 
sense the elasticity of such materials is expected to be 
very different from classical rubbers and strong devia- 
tions from the non-Hookian deformation behavior 
(as observed in the low deformation regime in conven- 
tional rubbers) must be expected. Actually, the elastic- 
ity of such Olympic gels has not been calculated yet and 
we attempt to present a simple argument in this issue. 
Apart from the case of synthetic materials these consid- 
erations are important for biological systems too. It is 
well known that highly entangled DNA rings exists and 
play an important role in biology 

Unless like in conventional soft materials, such as rub- 
bers the precise calculation of the modulus of topological 
gels is difficult because the exact topological state of the 
gel must be known. Specifying topological states is a 
general problem in polymer physics and has been first 
discussed in The corresponding mathematical prob- 
lem is the classification of knots and links |^ . Despite its 
incompleteness already noted in Q , in general the Gauss 
invariant is used in polymer physics because it is the most 
simple invariant that explicitly contains the polymer con- 
formation in contrast to algebraic invariants in knot the- 
ory. Even using the Gauss invariant, it is, in general, 
not very simple to describe the linking status of the net- 
work. Only for the easy case of a non-concatenated melt 
of rings this seems to be possible, but already here many 
complications appear |10[| . 

The synthesis and preparation of such gels is also a 
very difficult task. Some of the problems are already 
discussed in de Gennes' book [Q. The crucial point is 
that the ring closure of the chains must be carried out at 
concentrations larger than the overlap concentration c*. 
The reaction process is carried out in two steps: First a 
certain amount of rings is cyclicized by end group reac- 



tion. Then further linear chains are added to the given 
sample. These additional rings are then cyclicized again. 
With variation of the concentration different entangle- 
ment numbers can be expected. The solvent must then 
be evaporated to get the bulk network. In such an ideal 
synthesis the topological state of the network depends 
on the concentration of preparation. One limit is the cy- 
clizaton in the linear melt state (a gedanken experiment). 
Such prepared topological gels will be important for moti- 
vation of the following scaling analysis. Assume therefore 
a condensed melt of linear chains. It is well known, that 
the excluded volume forces are screened out [THTl|]. The 
screening of the excluded volume can also be asserted to 
the strong interpenetration of linear (= one dimensional 
connectivity) chains As a result, the chains behave 
gaussian in the melt. If then each long chain is closed to 
a ring [in Gedanken), the size of the ring polymer would 
not change. The reason for this is, that each ring is (in 
average) concatenated by n ~ q6'N^^'^, where q is the 
density of the melt, i the typical size of a monomer, and 
N the degree of polymerization. In the limit of scaling, 
n can be identified by the winding number. Such natural 
conjectures have been also mentioned by Gates et al. 
For a more rigorous definition of the winding number in 
such gel-like systems we refer the reader to our previous 
publication [ p^ . In this paper we had already confirmed 
the scaling conjecture presented first by Gates et al. |l^] 
by a careful treatment of the Gaussian linking number. 
The main result there was to show that the conformation 
of a ring in an unlinked melt of rings is naturally non - 
gaussian, i.e., the size R oc iV", where v < 1/2. 

In the present note we treat dense linked Olympic 
gels without solvent. The discussion of solvent free gels 
yields the "bare" scaling of the elastic modulus without 
additional contributions from the excluded volume and 
swelling. 

To do this, we have to use non-conventional arguments, 
because most of the theories for ideal, ordinary rubber 
networks are "single chain theories", i.e. all the elas- 
ticity of the network is computed by the contribution 
of the elasticity of one single network strand, which is 
then multiplied by the number of elastically active chains 
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Q . These well known results can be summarized by the 
free energy of the network as function of deformation 
F cx Mk-oTX^ , where M is the number of active chains 
and A is the deformation ratio of the chain. Although this 
seems to be very simplified, the results are in reasonable 
agreement with low deformation experiments and recent 
simulations as long as higher order effects at large 
deformations are ignored. In this note we put forward a 
scaling theory in a similar spirit. To do this we estimate 
the single ring behavior and conclude from this basis the 
elasticity of the entire gel. For the subsequent analysis 
we need the following assumptions 

• The rings are not self-knotted 

• The topological state of the network can be de- 
scribed by a global winding number n, which on 
the average is the same for each ring in the net- 
work 

• The direct influence of excluded volume effects are 
ignored 

Again, these assumptions seem to be crude, but these 
turn out to be sufhcient for the type of analysis presented 
here. A mathematical formulation of the problem of elas- 
ticity of the Olympic gel will be presented elsewhere . 

In the following we start from a Flory type estimate 
of the free energy of an entangled ring in a Olympic gel, 
where the average winding number of each ring is as- 
sumed to be n. The free energy of an entangled ring in 
a network can be written as 



N i?3 
F = kBT—r + const- 
it^ nN 



(1) 



Here the first term is the gaussian elastic part of the ring 
in the gel and the second term represents the pressure ex- 
perienced by the ring that comes from all the surround- 
ing rings. If the ring is not entangled with the others the 
second term would read B? /N obtained by replacing n 
by n -I- 1 in Eq.(|l|) and n = 0. Note that in this case 
the present ansatz for the free energy agrees with the 
one proposed by Gates and Deutsch for non concate- 
nated ring melts. If, however, other polymers are entan- 
gled with the ring under consideration they exert a re- 
pulsive force per winding number from inside and outside 
the ring, and consequently they do not contribute to the 
packing pressure. By this, we mean that the dense melt of 
rings is as closely packed as possible. Thus one ring expe- 
riences a pressure induced by the surrounding ones. The 
packing term is, however, reduced by the average wind- 
ing number. The factor 1/n in front of the second term 
has its reason in the "screening" of the packing pressure. 
Each other ring, that is entangled with the ring polymer 
under consideration reduces the packing pressure. In the 
limit of the cyclization in the melt the free energy con- 
tribution must be of the order of 0(1), since the rings 
must be gaussian. Note that the free energy described 
by Eq.(p|) does not contain an upper critical dimension. 



The interaction term is therefore important in all spatial 
dimensions. Instead of an upper critical dimension it con- 
tains a limiting winding number such that R cx VlV, cor- 
responding to the melt cyclization (^edan/cenexperiment. 
Minimization of the free energy with respect to R yields 
the size for the ring in the network 



Rn cx n^/^iV^/s 



(2) 



where all (to the purpose of this paper) irrelevant con- 
stants have been dropped. The latter equation has inter- 
esting consequences. For small n the scaling results agree 
with those proposed by Gates et al. The rings ap- 

pear compressed in the melt. Similar findings have been 
put forward by a more rigorous theory by us recently 
I p^ and by numerical simulations of the problem ||l6|,|l^ . 
On the other hand, if the mean winding number is close 
to that given by the average density of the system, i.e., 
n cx VlV, the ring is gaussian at melt cyclization con- 
ditions. In this case the topological effects are screened 
and the ring finds itself in a natural melt environment. 
If the winding number is larger, for example n (x N than 
too many rings are connected with each other and the 
ring configuration stretches out. So far the requirements 
have been satisfied by the simple scaling ansatz. 

The next task is now to find the asymptotic distri- 
bution function for the typical size r of the ring in the 
gel. A reasonable form of the distribution function is 

given by P{r,N,n) — Mr^ f (^-^^, where A/" is a suit- 
able but uninteresting normalization. 6 determines the 
short distance behavior and the function f{x) has to be 
determined. In fact, for the elasticity the knowledge of 
the scaling function / is sufficient, because the power in 
front of the the distribution function yields only irrele- 
vant (logarithmic) corrections. To find the appropriate 
distribution function, we start from the asymptotic be- 
havior of the scattering function defined by the size of 
the ring. Then, by standard methods, such as steepest 
descent Fourier inversion 0| it is easily found that the 
asymptotic form of the distribution is given by 



P(r, n) cx exp < — 



r 
Rn 



5/3" 



(3) 



This asymptotic form will be sufficient in the scaling 
limit. This is indeed the key equation of the paper. To- 
gether with the equation ^ it yields the correct asymp- 
totic behavior of the distribution function. For low values 
of the average winding number n it contains the limits 
suggested in jl^] 

The Olympic gel is considered to be entropy elastic. 
Thus we may conclude that the elastic free energy of the 
typical ring in the Olympic gel as a function of its elon- 
gation r is given by 



F,{r)^T 



^5/3 



(4) 
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The tension is given by the derivative and thus we find 
a non Unear force extension relationship due to the non 
gaussian structure of the rings in the network, i.e.. 

The latter result is the force extension law for one ring in 
the Olympic gel, and corresponds to the according single 
chain deformation law in conventional rubber theory, i.e., 
/ = {T/N)r Therefore (non-)Hookian deformation 
behavior in Olympic gels cannot be expected. 

At this point we have the possibility to observe two 
different deformation processes and regimes. At low de- 
formations the entanglements do not act as severe con- 
straints, but have many degrees of freedom, similar to en- 
tanglements in conventional rubbers. It has been shown 
that there the entanglements yield a "softening" of the 
modulus, if compared to the classical gaussian theory 
P, p^pS| ] . This softening of the modulus corresponds to 
the slippage and sliding of entanglements. This happens 
if the deformation of the individual ring in the gel is such 
that only the mean contour between two entanglements, 
i.e., N/n take part on the deformation. This defines 
the deformation ratio Aq = [N/n)/{RN) = {N/n^fl^. 
Note that Aq in the melt preparation conditions, when 
n oc -s/^V, is which is physically sensible: At 

high degree of entangling the topological constraints act 
immediately as crosslinks from the lowest deformation. 
From the arguments presented in ref. |§,|8| the maximum 
deformation can be estimated to Amax — 0{\/ N /n). The 
latter is always less than Aq as long as n < N^^'^ . For the 
validity of the scaling arguments this must always be the 
case here. 

Consider first the low deformation regime A < Aq. In 
the low deformation regime of Olympic gels the slippage 
of the topological constraints dominates. The main prob- 
lem with the low deformation regime is that the relevant 
chain length is not fixed |p|, p^|j2Q |. To see this point con- 
sider conventional rubbers, where the fixed length scale is 
given by the mesh size. In such Olympic gels a clear length 
scale (at least in the low deformation regime) cannot be 
defined, because the system is ruled by a large number 
of degrees of freedom. 

It has been shown that the deformation process can 
then be described by an effective distribution P{r) = 
/ dNp{N)P{r, N, n) function on the level of a single ring. 
As a consistent model we choose for p{N) the entangle- 
ments slack p{N) cx exp(— A^/A^o), which has been 
successfully applied to entanglement problems (see ||] 
and references therein). A'o is a mean excursion of the 
ring form the most probable conformation. The effective 
distribution is then given by the asymptotic form 



P(r, No^n) cx exp 



1/5^2/5 



(6) 



ensemble of Mr rings (per unit volume) is then given by 
multiplying the above equation by the number of con- 
straints present. These are the number of entanglements 
A/rti. To introduce the deformation we replace r by Ai?7v 
and average thus over the conformation as in the simplest 
theories in classical networks. The total free energy is 
then estimated by 



F = TM^nX 



(7) 



The non-Hookian linear increase of the free energy is en- 
tirely due to the large degrees of freedom of the con- 
straints, and is in some way similar to the low defor- 
mation regime in highly entangled rubbers, with strong 
entanglement sliding, when the sliplink contribution be- 
comes very weak (see ref. Q for details). Consequently, 
the free energy of the Olympic gel is linear in the elonga- 
tion A, and the interesting result is that the force needed 
for elongating the gel is constant, i.e., / = TtiMr in the 
low deformation regime. In terms of macroscopic vari- 
ables it is given by / = {Tq)/N, where q = M^N is 
the macroscopic density of the gel. Such a deformation 
regime is not observable in classical (highly entangled) 
gels, because in Olympic gels the degrees of freedom of 
the non crosslinked polymers are much larger, leading to 
a very weak solid. This is, of course, because Olympic 
gels consist only of entanglements. 

At larger deformations (A > Aq) the individual chains 
are deformed also, and the topological constraints act as 
crosslinks. Therefore the free energy must be propor- 
tional to the effective umber of constraints, i.e., Mrti. In 
this case we obtain 



F ^ MrtiA^/^ 



g-A5/3 



(8) 



which is consistent with eq. (^ since the mean size of the 
ring is not altered. The macroscopic free energy of an 



Obviously il/Rn is the effective number of crosslinks in 
the (afhne) deformation regime, where the topological 
constraints act almost as crosslinks. The measured force 
/ is then give by the derivative of the free energy with re- 
spect to the deformation A., i.e., / = Q{n/N)\^^^ and is 
larger by a factor of n due to the number of constraints 
but much weaker in the deformation dependence when 
compared to classical rubbers. In the latter case the 
force is roughly given by / = g{l/N)X The factor 
of l/{N/n) in the force can be indeed interpreted by the 
effective meshsize of the Olympic gel. distance 

In summary, we have presented a simple calculation 
of the elasticity of Olympic gels as a classical example of 
weak solids. We found two relevant deformation regimes 
which are determined by the topological state of the net- 
work. The first is the non-affine regime where the modu- 
lus is very weak indeed and the scaling is determined by 
the average winding number and the degree of polymer- 
ization of the rings. In the second deformation regime the 
topological constraints act similar as crosslinks. There- 
fore the (low deformation A « 1) modulus is given by 
Gaffinc — oin/N). When the mean winding number n is 
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of the order of \fN then the modulus becomes larger com- 
pared to the classical rubber. This case corresponds to 
ring closure of the polymers in the melt state. This high 
modulus in the affine regime is naturally determined by 



the large number of effective crosslinks. We expect that 
the present results have some applications in biological 
systems, too. 
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